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A theory of electrical transport through molecular wires is used to estimate the electronic fac- 
tor in the intramolecular electron transfer (ET) in porphyrin-nitrobenzene supermolecules, and 
to analyze its structure. The chosen molecules have complex donor and acceptor configurations, 
and relatively simple structure of the bridge, which enables us to concentrate our studies on 
donor/acceptor coupling to the bridge. We present analytical and numerical results concerning 
the effect of donor/acceptor coupling to the bridge on the ET process in molecules with complex 
donor/acceptor subsystems. 



(N 



o 

CO 



T3 

C 

O 

Q : 

(N 
> 

00 
(N 
^1- 
O 
(N 
O 
-i— > 

a 



i 

-a 
c 

o 
o 



X 



PACS numbers: 05.60.Gg, 36.20.-r 



Long-distance electron transfer (ET) plays a central 
role in many biological processes [lj. It has been es- 
tablished that molecular ET is essentially a combination 
of nuclear environment fluctuations and electron tunnel- 
ing 2]. The expression for the ET rate including both 
electronic and nuclear factors was first proposed by Mar- 
cus 0, 0, @ • The existing theories of long-range ET are 
mainly concentrated on the analysis of the electron tun- 
neling through the bridge (Refs |S II H IS ES IH E3 ) • 
Hitherto, less attention has been paid to systematic stud- 
ies of donor/acceptor coupling to the bridge and its con- 
tribution to the electronic transmission coefficient. Usu- 
ally, theoretical models of ET present both donor and 
acceptor as single sites with one state per site. In realis- 
tic biological systems, however, donor and acceptor sub- 
systems are complex and include numerous sites which 
provide effective coupling to the bridge. 

In this letter, we study the effect of the donor/acceptor 
subsystems on the characteristics of the long-range ET 
in macromolecules. We use an apparent close simi- 
larity between this process and an electron transport 
through molecular wire placed between metallic con- 
tacts. This enables us to apply the methods developed 
m Refs. |H H 03 to calculate the electron trans- 
mission function T(E) which we can further analyze to 
separate out contributions from different donor/acceptor 
sites participating in the ET process. We also calculate 
current- voltage characteristics of the intramolecular elec- 
tron transport. 

In further studies, we simulate our donor and accep- 
tor subsystems as sets of semi-infinite tight-binding ho- 
mogeneous chains. Each chain is attached to a site of 
donor/acceptor which can be effectively coupled to the 
bridge. To simplify the following calculations we assume 
that the chains do not interact with each other. This 
model for the donor/acceptor subsystem is a general- 
ization of that adopted in the th eory of electron trans- 
port through molecular wires [13, [lj, [lJj, [lg where each 
metallic electrode was treated as a single chain. How- 
ever, in the present analysis we need this generalization 
to take into account the complexity of the donor/acceptor 



subsystems. 

To proceed we consider each element of the donor- 
bridge-acceptor system as an array of potential wells. 
Each well represents one atomic orbital of the cor- 
responding part of the molecule (donor, acceptor or 
bridge). These wells are sites of the donor-bridge- 
acceptor system. Within this approach, any atom will be 
represented by a set of sites corresponding to its states, 
and we do not distinguish between sites corresponding 
to the same atom, and those corresponding to different 
atoms of our molecular system 0, 0] . 

In the following calculations we start from a tight- 
binding Hamiltonian for the bridge: 



H = Hn 



Hi 



(1) 



whose matrix elements between states \i > and \ j > cor- 
responding to the i-th and j-th sites are given by: 



(2) 



where 5ij is Kronecker symbol, Vy = when i = j and 
only states associated with valence electrons are consid- 
ered. The diagonal matrix elements evi are ionization en- 
ergies of electrons at sites i, while the off-diagonal matrix 
elements Vy = Vji are the hopping strengths between the 
ith and jth sites. Both direct and exchange energy con- 
tribution are included in Vy. As well as it was carried 
out in 0, ^| we also take into account self-energy cor- 
rections arising due to the coupling of the donor (Ep) 
and acceptor (E^) to the bridge. As a result, we arrive 
at the effective Hamiltonian for the bridge, i.e., 



Heff = H + H 



V 



(3) 



Generalizing the results of Ref. 18] for the case when 
we have several donor and acceptor sites coupling to the 
bridge, we obtain 



T>1 



E - e k - o-fc 



{H A ) l0 = (£, 



A 2 , 

-ei — cTi 



(4) 



(5) 
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Here, T>k% and An are, respectively, coupling strengths 
between the fc-th donor site or the Z-th acceptor 
site and the i-th site of the bridge, and ak,l = 

\ \^k,i — ^\J^lk,i ~ ®k 2} are tne self-energy corrections 
of the semiinfinite chains attached to the corresponding 
sites The parameters Ok.i = E — ek,i, tk,i and jk,i are 
the ionization energies of electrons at the corresponding 
donor/acceptor sites, and the nearest-neighbor hopping 
integrals for the chains. Summation in Eqs. (@J and (JSJ) 
is carried out over all donor/acceptor sites coupled to the 
bridge. Due to the presence of the self-energy corrections 
^t>,A: the eigenvalues of the effective Hamiltonian in- 
clude imaginary parts I\ which represent broadening of 
the bridge energy levels Ei originating from the coupling 
of the bridge to the donor and acceptor systems. The 
energy levels are broadened further if include scattering 
processes in the bridge. Here, however, for simplicity, we 
separate the effect of electronic dephasing, assuming that 
coherent electronic tunneling through a bridge to be the 
predominant mechanism of ET. 

Treating the long-range ET processes, we can neglect 
the direct coupling of donor sites to acceptor sites. Then 
the probability amplitude for the transition between the 
fc-th donor site and the l-th. acceptor site is given by 0] 



Tu — 7 ^kiGijAji 



(0) 



Here, the sum is carried out over the bridge sites, Gy is 
the matrix element of the Green's function corresponding 
to H eff , i.e., 



Gij = 



(e - H eff y 



(7) 



An electric tunneling current / flowing from donor to 
acceptor through the bridge in the presence of a small 
applied voltage V has the form [2(| 



I : 



2ire 



k,l 



f (E k ) [1 - / (Ei + eV)] I Tl, I 6 (E k - E t ) 



(8) 

where f(E) is the Fermi function; e is the charge of 
the electron. Assuming that the applied potential varies 
linearly in the molecule, we can use the approximation 
V = Vo(Lb/ L mo i), where L}, and L mo i are the lengths 
of the bridge and the whole molecule, respectively, and 
Vb is a voltage applied across the whole molecule. Start- 
ing from the expression (8) and following a usual way we 
arrive at the standard formula fell: 



e r°° 

7 = — / dE T(E)[f (E - Ml ) -f(E- M2 )] 



(9) 



where the chemical potentials fj,± and /12 are determined 
by the equilibrium Fermi energy of the bridge Ep and 
the effective voltage V across the bridge 0]j i.e., 



Mi 



E F + (1- r?) eV; i i 2 = E F ~ r/eV. 
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FIG. 1: The calculated current (nA) - voltage (V) charac- 
teristics for different temperatures. 



Here the parameter rj characterizes how the voltage V is 
divided between the two ends of the bridge. The equi- 
librium Fermi energy for the bridge can be determined 
from the equation (See Ref. (lfj): 



N 



arctan • 



E f - Ei 



r, 



(10) 



where N is the number of electrons on the bridge part of 
the molecule. 

The electron transmission function included in the ex- 
pression (9) is given by the formula 



T(75) = 2^A i |G i ,fA i 



(11) 



where the summation is carried out over the states of the 
bridge, and the quantities A,.j arc defined as 



A i =Im(S 2? ) 



This formula agrees with the result obtained earlier in 
the study of electron transport through molecular wires 

We see that contributions to the transmission func- 
tion T(E) from different donor and acceptor sites are 
determined by the factors A,-j. ft follows from Eqs. (@J 
and JSJl that these factors strongly depend on energy- 
dependent denominators. When the tunneling energy E 
takes on a value close to e/^j, the corresponding term 
in Eqs. and (JSJ can surpass all remaining terms. 
Thus, for different values of E, different donor and ac- 
ceptor sites can predominate in the ET process. Ex- 
pressions (0J and JSJ) enable us to determine these dom- 
inating sites for any interval of tunneling energy E. In 
Fig.l we present results for the calculated I-V charac- 
teristics for the porphyrinc-nitrobenzene macromolecules 
for different temperatures (T = IK, 50K, 100K, 300K). 
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FIG. 2: The calculated electron transmission function 
through the bridge T(E) within a certain range of energies 
E(eV). 



In calculations we estimate Ep using the expression (10) 
and we choose T) = 0.3. Matrix elements of the bridge 
Hamiltonian (1) as well as the ionization energies for 
the donor/acceptor sites are obtained from the extended 
Huckel method. The parameters 7^; representing hop- 
ping integrals for donor/acceptor chains are estimated as 
averages of matrix elements which correspond to the cou- 
pling of these sites to their nearest neighbors within the 
donor/acceptor part of the molecule. 

Comparison of our results with those obtained for I- 
V curves characterizing the electric transport through 
molecular wires [lij shows their resemblance at room 
temperatures. At low temperatures we notice that our 
I-V curves exhibit a step-like behavior and the plateaus 
get narrower as the temperature raises. Associating the 
low-temperature I-V curves with the result for the elec- 
tron transmission function T(E) shown in part in Fig. 2 
we see that the steps reflect the structure of T(E). The 
latter is a set of peaks separated by the intervals where 
the electron transmission function takes on values close 
to zero. These peaks arise due to the contributions from 
different donor-acceptor sites as well as from the Greens 
function corresponding to the effective Hamiltonian for 
the bridge. At a given voltage V the difference of the 
Fermi functions in expression (9) differs from zero only 
in the interior of a certain energy range including Ep. 
Therefore, the tunneling current I is fully determined by 
the contributions of peaks of the transmission function lo- 
cated there. When we increase the voltage applied across 
the bridge the width of the relevant energy range en- 
hances, and extra peaks of the function T(E) contribute 
to the current. When a new peak appears this gives rise 
to a sudden enhancement of /, so that the step exhib- 
ited in the I — V curves reflect the structure of T(E). 
At higher temperatures Fermi distribution functions in 



the integrand of (9) lose their step-like character and the 
plateaus are washed out. Another reason for peaks in 
the electron transmission function (and the steps in the 
I-V curves, respectively) to be eroded is the electronic 
phase-breaking effect which arises in complex molecules 
due to stochastic fluctuations of ion potential. We con- 
jecture, however, that at low temperatures the electronic 
dephasing reduces to same extent, so that the structure 
of T(E) can be revealed. 

In summary, we have presented a theoretical model 
for the intramolecular ET which resembles those used 
to study the electric transport through molecular 
wires. Our model enables us to analyze in detail the 
donor/acceptor coupling to the molecular bridge which 
is important for further investigations of long-range ET 
through macromolecules. 
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